Nucleation mechanism for the direct graphite-to-diamond phase transition 
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Graphite and diamond have comparable free 
energies, yet forming diamond from graphite is 
far from easy. In the absence of a catalyst, pres- 
sures that are significantly higher than the equi- 
librium coexistence pressures are required to in- 
duce the graphite-to-diamond transition 1 7 . Fur- 
thermore, the formation of the metastable hexag- 
onal polymorph of diamond instead of the more 
stable cubic diamond is favored at lower temper- 
atures^—. The concerted mechanism suggested 
in previous theoretical studies^ — cannot explain 
these phenomena. Using an ab initio quality 
neural- network potential 13 we performed a large- 
scale study of the graphite-to-diamond transition 
assuming that it occurs via nucleation. The nucle- 
ation mechanism accounts for the observed phe- 
nomenology and reveals its microscopic origins. 
We demonstrated that the large lattice distor- 
tions that accompany the formation of the dia- 
mond nuclei inhibit the phase transition at low 
pressure and direct it towards the hexagonal di- 
amond phase at higher pressure. The nucleation 
mechanism proposed in this work is an important 
step towards a better understanding of structural 
transformations in a wide range of complex sys- 
tems such as amorphous carbon and carbon nano- 
materials. 

Static compression of hexagonal graphite (HG) results 
in the formation of metastable hexagonal diamond (HD) 
at temperatures around 1200-1700 K— £~— and cubic di- 
amond (CD) at higher temperatures^ 3 - - — i Although 
the transition pressure is sensitive to the nature of the 
graphite samples neither of the diamond phases has been 
observed to form below ^12 GPa. This pressure is 
significantly higher than the graphite-diamond coexis- 
tence pressure approximated by the Berman-Simon line 
P(GPa) - 0.76 + 2.78 x 10~ 3 r(K)H. 

Despite being an area of intense theoretical re- 
search^— the microscopic mechanism of the formation 
of metastable HD and the reason for the remarkable sta- 
bility of graphite above the coexistence pressure are still 
unknown. Computer simulations, which could help re- 
solve these issues, have been hindered because of the in- 
ability of empirical potentials to describe the energetics 
of the transformation accuratel y 13 ' 15 and the computa- 
tional expense of more reliable ab initio methods. In the 
latter case, short simulation time and small system size 



(i.e. several hundred atoms) force the transition to occur 
in a concerted manner with the ultrafast (~ 10~ 2 — 1 ps) 
synchronous formation of all new chemical bonds accross 
the entire simulation bo x 11 ' 12 ' 16 . While concerted mech- 
anisms can be observed at shock compression^—, the 
transformation under static conditions is expected to pro- 
ceed via nucleation and growth. 

It has been estimated that because diamond has an 
extremely high surface energj*^ its critical nuclei may 
contain thousands of atoms 20-22 . Hence, tens or even 
hundreds of thousands of atoms are required for model- 
ing the diamond nuclei and the surrounding graphite ma- 
trix. Direct ab initio simulations of systems of this size 
are outright impossible. Therefore, theoretical studies 
of the nucleation have been restricted to simple contin- 
uum models 2 ^— , which neglect the anisotropic nature of 
graphite, use significantly different estimates of the sur- 
face energy terms, and ignore the distortion of graphite 
around the growing nuclei. 

In recent works, we have demonstrated that high- 
dimensional neural networks (NN)^ 3 - are capable of creat- 
ing accurate representations of ab initio potential energy 
surfaces of numerous element s 13 ' 23 ' 24 . Even in the case of 
graphite and diamond which are very differently bonded, 
a NN potential predicts all relevant properties in quan- 
titative agreement with ab initio and experimental data 
(see Methods and Ref . [Ill) . Furthermore, the computa- 
tional efficiency of MD based on the NN potential enables 
us to extend time and length scales accessible to simu- 
lations and, thus, to perform the first atomistic study of 
homogeneous diamond nucleation from graphite. 

The energetics of the nucleation was studied at zero 
temperature by seeding diamond nuclei of various sizes 
inside a periodic ^100x100x100 A graphite matrix con- 
taining ^145,000 atoms (see Methods). Hexagonal and 
rhombohcdral graphite (RG) lattices were used as ini- 
tial structures for the formation of HD and CD nuclei, 
respectively, because of high symmetry of the HG— >-HD 
and RG— i-CD transformation pathway s 8 ' 10 . HD and CD 
nuclei were formed by shifting atoms in the graphite lat- 
tice in the manner shown in Fig. [T]A_. Buckling of basal 
planes into the "chair" conformation in the RG lattice 
leads to CD whereas the "boat" buckling of the HG lat- 
tice results in the formation of HD. Such distortions of 
graphite planes have been observed in previous ab ini- 
tio simulations 1 ^ and satisfy the orientation relations be- 
tween graphite and diamond crystals discovered experi- 
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FIG. 1. A. Pathways for RG— s-CD (green) and HG— ^HD (red) transformations. OG denotes an (unstable) intermediate 
orthorhombic graphite phase. B. DFT (lines) and NN (circles) enthalpy profiles for the concerted mechanism: RG— >CD 
(green) and OG— >HD (red). The enthalpy of undistorted RG is taken as zero. 
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(100) HD || (001) G and [100} HD || [100] G , 
(1U) CD || (001) G and [U0] CD || [100] G . 

Our calculations show that the nucleation process is 
strongly influenced by pressure. Both the nucleation bar- 
rier and the size of the critical nuclei decrease rapidly as 
the pressure is increased (Fig. [5] and [3]) . Examination of 
the atomic structure of the diamond nuclei reveals the 
microscopic origin of this phenomenon. Diamond nuclei 
are generally highly non-spherical and contain a diamond 
core (red atoms in Fig. [5]) separated from a relaxation re- 
gion of graphite lattice by a thin (~5 A) high-enthalpy 
interface. To form interlayer bonds in the core, large 
regions in a number of graphite layers have to be buck- 
led and bent in the [001] G direction. At low pressure 
these distortions are energetically costly because of the 
significant mismatch between the [001] G lattice parame- 
ters of the parent phase and the nuclei. At 10 GPa, we 
did not observe the formation of interlayer bonds even 
around large diamond seeds. At higher pressure the mis- 
match decreases and the diamond cores become increas- 
ingly more spherical (Fig. [5]). The same mismatch is the 
reason for the [001] G distortion of the graphite crystal 
around the nuclei. The extent of the distorted region 
decreases with the increase in pressure. 

Since the size of the distorted region around the dia- 
mond core grows linearly with the size of the core the 
free energy of the diamond nucleus in graphite can be 
written as2&22 

AG = (Ag s + Ag^V + <rS, (1) 

where Ag^ is the difference between the free-energy den- 
sities of bulk diamond and graphite, Ag s is the positive 
misfit strain free-energy per nucleus volume and a is the 
interfacial free-energy density. S and V are the surface 



area and volume of the nucleus, respectively. Due to 
the considerable mismatch between the lattice parame- 
ters, the large strain energy term outweighs the relatively 
small bulk term at P < 10 GPa and makes nucleation im- 
possible (i.e. the volume term becomes positive). As a 
result the 10 GPa nucleation curve does not show any 
signs of flattening out even at R — 50 A. In contrast, 
the 20 GPa curve is predicted to reach the maximum 
around 25-31 A and 560-630 eV. Hence, our results are 
consistent with the observed ^12 GPa minimum pres- 
sure threshold for diamond formation. The negative bulk 
term Ag^ increases in magnitude with pressure (Fig.[Tj3) 
while the positive strain term becomes smaller (i.e. note 
the extent of the relaxation region around the diamond 
core in Fig. [3]) leading to the observed decrease in the 
nucleation barriers (Fig. 13]). 

Comparison of the transition barriers for the con- 
certed pathways and the nucleation process (Table QJ 
shows that the nucleation is energetically more favorable 
in the pressure range used in compression experiments. 
Hence, nucleation represents a more realistic mechanism 
for diamond formation. However, at higher pressure 
(~ 50 GPa) the graphite crystal approaches the lattice 
instability point and the activation barriers for the con- 
tinuous transformations are lower (Fig.Qj3 and Table HJ. 
In this high-pressure limit, diamond domains can appear 
spontaneously throughout the graphite matrix without 
the formation of a well-defined graphite-diamond inter- 
face. For example, at 60-70 GPa, the formation of only 
several interlayer bonds in close proximity to each other 
is enough to initiate a rapid irreversible growth of a dia- 
mond crystal. 

The activation barriers for the concerted mecha- 
nism are determined by the energetics of buckling of 
graphite basal planes into the "chair" or "boat" configu- 
rations. The considerably lower activation barrier along 
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FIG. 2. Pressure dependence of the shape and size of the CD nuclei. The optimized critical nuclei are shown for 30, 40, and 
50 GPa whereas the 20 GPa nucleus is below the critical size and is presented to compare its shape with that of the 30 GPa 
nucleus. Atoms are colored according to the values of the tetrahedral atomic order parameter defined to distinguish graphite 
(Xi ~ 0) and diamond (x,i ~ 1) configurations (see Methods). 
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FIG. 3. Pressure dependence of the barriers for the RG— >CD 
(circles) and HG— >HD (squares) diamond nucleation. The 
20 GPa curve is predicted to reach the maximum around 25- 
31 A and 560-630 eV. The sharp drop of the curves at high 
pressure provides evidence that there are extremely low barri- 
ers for adding atoms to the critical nuclei that prevent stabi- 
lization of supercritical CD nuclei even at low temperatures. 



the "chair" buckling mode favors the formation of CD 
(Fig.QjB and Table[TJ while in experiments the metastablc 
HD phase is often observed. This disagreement between 
the theoretical predictions and experiment is resolved if 
the nucleation mechanism of diamond formation is con- 
sidered. Our calculations show that the nucleation barri- 
ers for the HG— >HD and RG— >CD processes are not very 



TABLE I. Comparison of the enthalpy barriers of the con- 
certed pathways and nucleation. 



Pressure 
(GPa) 


Nucleation^ 
(meV/atom) 


Concerted (meV/atom) 
CD b HD b 


30 


70-90 


130 185 


40 


40-60 


80 140 


50 


110-280 


50 93 



a Calculated by dividing the nucleation barrier by the number of 
atoms in the diamond core (xi > 0.8) of the critical nucleus. 
This value gives the upper bound on the nucleation barrier. 

b Pathways are shown in Fig. [TK with green (CD) and red (HD) 
colors. 



different (Fig. |3|). The similar heights of activation bar- 
riers is the consequence of close energy of CD and HD 
nuclei and their surface energies^. 

In the nucleation process, the graphite basal planes 
undergo local "in-layer" distortions that bring atoms 
into appropriate stacking sequence inside diamond nu- 
clei. The slightly higher barriers along the pathway lead- 
ing to the HD phase are a consequence of this local [210] G 
distortion. This distortion brings the neighboring layers 
from the hexagonal (AB) to orthorhombic (AB') stack- 
ing (Fig. [OA.) inside the HD nuclei. The HG^CD nu- 
cleation process is considerably more difficult than both 
the RG— ?>CD and HG->HD nucleation. The strains nec- 
essary to distort the lattice from the AB stacking in HG 
to ABC stacking inside CD nuclei prevent stabilization 
of the CD nuclei inside the HG lattice in our simulations. 
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In-layer distortions have not been observed in any 
of the previous theoretical studies of this phase transi- 
tio n 11 ^ 12 ' 16 because in-layer stresses result in the artificial 
sliding of graphite layers in small-cell simulations. Our 
results suggest that the in-layer distortions in large-cell 
simulations determine the energetics of the transforma- 
tion and the overall direction for the phase transition. 
The higher in-layer strains for the HG— >HD transforma- 
tion relative to the HG— >-CD transition explain why HD 
is formed from graphite, which mostly consists of the 
hexagonal form. On the other hand, our studies show 
that the formation of CD is favored in the RG lattice be- 
cause no in-layer distortions are required for the RG— >CD 
process. This implies that CD nucleates in the HG sam- 
ples around stacking faults (e.g. [AB] n [C4] m ), which 
are common defects in the HG lattice. High temperature 
fluctuations can also activate the sliding of graphite basal 
planes in HG with the formation of the RG regions and, 
hence, favor the nucleation of the CD phase. 

In conclusion, our study of the homogeneous diamond 
nucleation offers new insights into the atomistic mecha- 
nism for the direct graphite-to-diamond phase transition. 
We have demonstrated that the transformation does not 
occur below ^10 GPa in the static compression experi- 
ments because of the prohibitively large strains accom- 
panying the formation of diamond nuclei at low pressure. 
We have also shown that, at higher pressures, the nucle- 
ation mechanism is favored over the concerted transfor- 
mation. Larger distortions of the graphite lattice around 
the CD nuclei compared to those around the HD nu- 
clei explain the formation of the metastable HD phase. 
At yet higher pressures, the transition is continuous and 
proceeds without formation of a well-defined graphite- 
diamond interface. 

A finite temperature MD study of the thermodynam- 
ics and kinetics for heterogeneous nucleation on defects in 
the graphite lattice will be carried out in the future. The 
computational and theoretical models presented here of- 
fer new opportunities for investigating of complex struc- 
tural transformations in a wide range of carbon-based 
materials. 

Methods. The NN potential was created to repro- 
duce the Perdew-Burke-Ernzerhof (PBE) density func- 
tional^ energies for both phases in the pressure range 
up to 100 GPa, as described in our previous paper—. 
The PBE functional in combination with a dispersion 
corrected atom centered pseudopotential 30 accounts for 
the van der Waals interactions in graphite and describes 
experimental structural, elastic and vibrational proper- 
ties of the diamond and graphite phases^. The NN fit- 
ting procedure introduces only a small error (RMSE of 
the independent test sets is 4.0 meV/atom) in addition 
to the numerical (convergence) error of the DFT calcula- 
tions. Thus, the NN-potential is expected to describe all 
properties and processes in graphite and diamond with 
an accuracy comparable with that of the PBE functional. 

Fig. [TJ3 shows that the energetics of concerted trans- 
formations in small cells obtained with the NN poten- 




C-C distance, d (A) 

FIG. 4. DFT (lines) and NN (circles) enthalpy profiles for a 
nucleation pathway of CD. The enthalpy of undistorted RG 
is taken as zero; the 20 and 40 GPa curves are shifted up for 
clarity of presentation. 

tial is indistinguishable from the density functional the- 
ory (DFT) results. In quantitative agreement with the 
ab initio calculations of Tateyama at al^, the stability 
of diamond relative to graphite increases with pressure 
whereas the barrier separating two phases decreases. At 
a pressure of 80 GPa graphite reaches a lattice insta- 
bility point (i.e. one of the barriers becomes zero) and 
undergoes an ultrafast transformation to diamond as was 
previously observed in ab initio simulations by Scandolo 
et aiii. The NN-driven MD simulations at constant pres- 
sure and temperature correctly reproduce the mechanism 
of this concerted transformation and are in perfect agree- 
ment with the results of Scandolo et al. 

Fig. [4] demonstrates that the formation of a diamond 
nucleus in a small system is also excellently described by 
the NN. 

The energetics of the nucleation of diamond from 
graphite was studied at zero temperature on a model 
system of 145,152 carbon atoms arranged in a graphite 
lattice in a periodic ^100x100x100 A simulation cell. 
Diamond nuclei of various shapes and sizes were seeded 
by shifting atoms within a certain predefined region of the 
cell in the directions shown in Fig. Q~|A. and constraining 
distances between appropriate pairs of atoms to the val- 
ues corresponding to the C-C distance in diamond. RG 
and HG was used as a starting lattice for the formation of 
CD and HD nuclei, respectively. 30 ps constant pressure 
MD simulations at 1000 K were performed to relax atoms 
around the constrained region followed by quenching to 
zero temperature. Finally, the constraints were removed 
and the geometry was optimized at constant pressure to 
obtain fully relaxed nuclei. We found that regardless of 
the shape of the initial constrained region all relaxed nu- 
clei have shapes similar to those shown in Fig. [3] We 
verified that all simulations results are converged with 
respect to the system size by performing additional test 
calculations for a ~200x 200x200 A simulation cell. 

The tetrahedral order parameter \i f° r atom i is de- 
fined as follows: Xi = n^-i) J2j>k nijn ik (\ + cos9j ik ) 2 , 
where the cutoff function is n,*j = [l + e( riJ_ro - | / A, ']~ 1 , the 
atom coordination number is ni = ^2 ■ riij and tq = 1.6 A, 
Ar = 0.005 A. 
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